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I . Abstract 

| We study the long time asymptotics of a nonlinear, nonlocal equation used in the modelling 

5_i ■ of granular media. We prove a uniform exponential convergence to equilibrium for degenerately 

convex and non convex interaction or confinement potentials, improving in particular results by 
J. A. Carrillo, R. J. McCann and C. Villani. The method is based on studying the dissipation 
OO ■ of the Wasserstein distance between a solution and the steady state. 
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Introduction 

We consider the problem of convergence to equilibrium for the nonlinear equation 

d m = Am + v-(m(vv + vw*m)) t>o,xeR n . (l) 



This equation preserves mass and positivity and we shall be concerned with solutions which 
are probability measures on M. n at all times. It is used in the modelling of space- homogeneous 
granular media (see [3]), where it governs the evolution of the velocity distribution Ht{x) of a 
\ particle under the effects of diffusion, a possible exterior potential V and a mean field interaction 

through the potential W; we shall keep the variable x instead of v (for the velocities) for 
notational convenience. 

Steady states may exist as a result of a balance between these three effects, and we are 
concerned with deriving rates of convergence of solutions towards them. Following [3J, this issue 
has raised much attention in the last years and has been tackled by a particle approximation 
and logarithmic Sobolev inequalities in [15], by an entropy dissipation method in 0[T2] and by 
contraction properties in Wasserstein distance in [8j [TTj [5] (see also [9J [6] for related works in 
dimension one). The entropy method is based on studying the time derivatives of a Lyapunov 
function F of the equation (called entropy or energy), on the interpretation due to F. Otto of ([!]) 
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as a gradient flow of F (see [HI [T]) and on the notion of convexity for F due to R. J. McCann 
(see H5]). 

When V and W are uniformly convex, solutions converge exponentially fast to equilibrium, 
but the case of interest of [3] is V = and W(x) = \x\ 3 , whose convexity degenerates at 0. For 
this case, only a polynomial rate, or exponential but depending on the initial data, was obtained 
in [3 El [TT] . In the present paper we prove a uniform exponential convergence in Wasserstein 
distance of all solutions to the steady state. The method, introduced in the linear case in [JJ, is 
based on comparing the Wasserstein distance with its dissipation along the evolution. 

In Section[T]we derive the dissipation of the Wasserstein distance between solutions and easily 
deduce the classical contraction results. Section [2] is devoted to cases when the convergence is 
driven by the interaction potential W, with or without exterior potential V: in particular we 
prove the first result of uniform exponential convergence to equilibrium for degenerately convex 
interaction potentials and no exterior potential. In Section we give conditions to get an 
exponential convergence to equilibrium with both potentials being non convex. 

1 Dissipation of the Wasserstein distance 

Let P2(W l ) be the set of Borel probability measures on W 1 with J K „ \x\ 2 d/j, < oo. The Wasser- 
stein distance between two measures [i and v in P2(M n ) is defined as 



where tt runs over the set of joint Borel probability measures on M. with marginals \i and v. 
It defines a distance on P2(W l ) which metrizes the narrow convergence, up to a condition on 
moments. In the present work, convergence estimates will be given in terms of this distance, 
but interpolation estimates can turn such weak convergence estimates into strong convergence 
estimates. By the Brenier Theorem, if \x is absolutely continuous with respect to the Lebesgue 
measure, then there exists a convex function <p such that V = u, that is, J^gdv = 
J*Rn <?(Vy) d/j, for every bounded function g; moreover 



and V<^*#z/ = /x for the Legendre transform ip* of (p if also v is absolutely continuous with 
respect to the Lebesgue measure. We refer to [101 [TJ [19] for instance for these notions. 

We shall assume that V and W are C 2 potentials on W 1 , respectively a and /3-convex with 
a, f3 Gl, in the sense that V 2 V(x) a and X7 2 W(x) ^ j3 for all x G R n , as quadratic forms 
on M. n . Moreover we assume that the interaction potential W is even and that both V and W 
satisfy the doubling condition 



for all x,y G M. n , and analogously for W. 

We shall consider solutions which are gradient flows in the space P2(M. n ) of the free energy 





V(x + y) < C(l + V(x) + V(y)) 



(2) 



Ffa) = [ 



log fidx+ Vdfj, + - W(x - y) dfj,(x) d/j,(y) 




(3) 
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as developed as follows in pQ: Let fio be an initial datum in P2(M. n ). Then, by |13} Ths. 4.20 
and 4.21] or [1, Th. 11.2.8], there exists a unique curve fi = (fi t )t G C([0, +oo[, P 2 (R n )), locally 
Lipschitz on ]0, +oo[, satisfying the evolution variational inequality 

1 d„ T ,, s . » „. . a + min{/3, 0) TT o . . 

--W$(jh,<t) < F{a) - F(jh) ^± W^a) 

for almost every t > and all probability measure a in the domain of F. For all t > the 
solution fi t has a density with respect to the Lebesgue measure. Moreover the curve fi satisfies 
the continuity equation 

dtlH + V-(lHVt) = 0, t>0,x£R n 

in the sense of distributions, where the velocity field vt satisfies 

-lHV t = Vfi t + fi t VV + fit (VW * fit). 

In other words ft is a solution to (pQ), and the curve fi = (fitjt will be called the solution with initial 
datum no e P 2 (R n ). Finally t ^ J \v t \ 2 dfi t £ L~ ([0,+oo[) so the curve fi :]0,+oo[->- P 2 (M n ) 
is absolutely continuous (see [TJ Th. 8.3.1]); moreover, if v is another such solution with initial 
datum uq and associated velocity field Wt, then by pjJJ Th. 23.9] or [1, Th. 8.4.7] 

7:4: W 2(fM,^t) = - I (yipt(x) - x) ■ v t (x) dfit(x) - I (V(ft(x) - x) ■ w t (x) du t (x) 

for almost every t > 0; here ipt is a convex function on M. n such that V^ptj^fit = v t an d 
Vipl^vt = fit- Then one can perform a "weak" integration by parts as in |14l Th. 1.5] or [HJ 
Lem. 13] to bound from above the right-hand side by 

(A<pt(x) + A(pt(y<pt(x)) -2n + (A(y(pt(x),vt) - Ax, fit)) • (V^t(x) - xj) dfi t {x). 



Here Aip is the trace of the Hessian of a convex map <p on M. n in the Alexandrov a.e. sense and 
A(x,fit) = VV(x) + VW * fit(x). Moreover, since VW is odd, the term involving W is 

(VW * v t (V(pt(x)) ~ VW * fi t (x)) ■ (Vipt(x) - x) dfi t (x) 
{VW(Vift(x) - Vipt(y)) - VW{x - y)) ■ (Vip t (x) - x) dfi t (x)dfi t (y) 

= 1 If (vw(y<p t (x) - v<p t (y)) - vw(x - y )) ■ (v<p t (x) - v^ t (y) -{x- y)) dft t (x)dfi t (y). 

2 J JM?n 

We summarize as follows: 
Proposition 1.1 ([lj) If (fit) t and (vt)t are two solutions to then for a.e. t>0, 

5~| w 2 2 G u *> I/ t) ^ -Jv,w(y t \in) 

where, for v = Vtpjfcfi (and Aip the trace of the Hessian of tp in the Alexandrov sense), 

J VjW (u\fi)=f (a<p(x) + A(p*(Vip(x)) - 2n + (W(V<p(x)) - VV(x)) ■ (Vip(x) - x)) dfi(x) 
+ 111 (VW(V<p(x) - V<p{y)) - VW(x - y)) ■ (V V (x) - V<p(y) - (x - y))dfi(x)dfi(y). (4) 
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For t > we can expect the solutions to have smooth densities, and to have equality in 
Proposition II, l\ but we shall be content with the inequality (see [H]). 

Considering the dissipation of the distance between two solutions provides simple alternative 
proofs of contraction properties in Wasserstein distance derived in [TT| . For that purpose we 
first notice that given /i and v absolutely continuous with respect to the Lebesgue measure, 
and V(p#n = then Aip + A(p*(y<p) — 2n ^ [i a.e. (see for example |141 Th. 1.5] and [U 
Lem. 2.5]). This inequality says that the diffusion part of the equation always contracts two 
solutions, as it is classical for the pure heat equation. Then: 

• Suppose that V and W are respectively a and /3-convex with a£l and /3 < 0. Then the 
term involving V in @ is bounded from below by a W^ip, v) and the term involving W by 



\V<p(x)-Vtp(y)-(x-y)\Uii(x)dti(y)=pW$(ii,v)-P j(Vi P -x)d f i > fiWlM 
since /3 < 0. Hence, for two solutions (nt)t and {yt)t of (HJ an d almost all t ^ 

IjWi&t, vt) < -(« + P) w%{m, ut). 

Then by the Gronwall lemma we recover the contraction property of [U Th. 5]: 

W 2 (JH, n) < e-^ 1 W 2 {fi , vq), t> 0. (5) 
Suppose that W is convex and that there exist p, C > such that for all e > 

(W(y) - VV(x)) • (y- x) > Ce p (\y-x\ 2 -e 2 ), x,yeR n . (6) 
Then, by the same argument, 

l±Wi(^u t ) < -^{2W 2 (^v t )-e 2 ). 
We optimize in e and integrate to recover the polynomial contraction of [8l Th. 6] 

W 2 (fH, vt) < (W^ p {^ v Q ) + ctY 1/P , t> 0. (7) 



• Suppose that V and W are respectively a and /3-convex with a E R and /3 ^ 0. Then, 
again by the same argument, the contraction result © holds for any two solutions with 
same center of mass, that is, such that in W 1 J^ n xdjit = J K n xdvt for all t ^ 0. This was 
also proved in [51 Th. 5]. 

• Suppose that V is convex and that ([6]) holds for W instead of V . If moreover the center of 
mass of each solution is conserved, that is, if J* Rn xdfxt = Ln xdfio for all i ^ (this is the 
case for instance if V = 0) then the polynomial contraction (J7j) holds for any two solutions 
with same (initial) center of mass, recovering [8j Th. 6] and Th. 4.1]. 

In the first case with a + (3 > the bound © ensures the existence of a unique stationary 
solution to (TT]) in P2(W l ), and the exponential convergence of all solutions to it. In the third 
case with a + j3 > 0, and if moreover the center of mass is preserved by the evolution, then for 
any m 6 R n this ensures the existence of a unique stationary solution to (TT]) in P2(M. n ) with 
center of mass m, and the exponential convergence to it of all solutions with (initial) center of 
mass m. 

The following two sections are devoted to the obtention of explicit exponential rates of 
convergence of solutions to ([1]) in non uniformly convex or even non convex cases, having in 
mind the degenerately convex potentials of [3] and the double well potentials of [18] . 
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2 Influence of the interaction potential 



In this section we study the case when W brings the convergence. 



2.1 No exterior potential 

We first assume that V = 0. Then the evolution preserves the center of mass, and a solution 
jit should converge to a stationary solution /Xqo only if the initial datum hq and fi^ have same 
center of mass, since 

xd/j, 00 (x)— / xdfio(x)= / xdfi OQ (x)— / xdfit(x) 

JK. n il" il" 

should converge to 0: for instance it is bounded by W2(fJ>t, fJ-oo)- We could also assume that 
V ^ 0, but that the center of mass is fixed by the evolution, which is all we use. But to simplify 
the statements we assume V = 0. 

When W is degenerately convex, with a pointwise degeneracy, for instance W(x) = |x| 2+£ 
with s > 0, then the contraction property holds only with polynomial decay rate, see the last 
example in Section [TJ Then in [7] the authors proved an exponential convergence to equilibrium, 
but not with a uniform decay rate, but rather depending on the free energy F of the initial 
datum. In this section we prove a uniform exponential convergence for such potentials. 

Theorem 2.1 Let V = and W a C 2 convex map on W n for which there exist R and K > 
such that 

V 2 W(x) ^ K if \x\ ^ R. 

Then for all m £ lR n there exists a unique stationary solution /i™ € P2(W l ) to ([T]) with center 
of mass m; moreover there exists a positive constant C such that all solutions (nt)t to ( f T | ). for 
an initial datum with center of mass m, converge to /i™ according to 

Proof 

< The existence of a stationary solution //^ G P2(W l ) with center of mass and a positive 
density satisfying /^(x) = Z~ 1 e~ w * fM ° c ^ is given by Proposition IA.lt proof of i., with any 
b < K/2; here Z is the normalizing constant. Then /i™ = — m) is a stationary solution 
with center of mass m. Now Proposition 12.21 and Remark 12.31 ii. below ensure the convergence 
estimate to /x™ since //™ = e~ u jZ with U = W * //™ convex and bounded from below. Unique- 
ness follows. l> 



Proposition 2.2 Let W be a C 2 convex map on M. n for which there exist R and K > such 
that 

V 2 W{x) > K if \x\ ^ R. 
Let fi G P2(W l ) have a continuous density e~ u for which there exists M such that 

sup sup {U(z) - U(x) - U(y)} < M. (8) 

\x-y\<2R z£[x,y] 

Then there exists an explicit positive constant C, depending only on K, R and M , such that 
for all measures v with same center of mass as [i. 
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Remark 2.3 Hypothesis ([8]) on U holds on any of the following two instances: 

i. U is C 1 and U(x) - 2R sup |VE%)| ^ —M for all x € W l 

\x-y\<R 

ii. U is C 2 , V 2 U(x) ^ a(x) with a(x) < and U(x) + 2R 2 inf a(y) ^ —M for all x; for 

\x—y\<R 

example, U is C 2 and bounded from below and V 2 U(x) a for all x and a constant a. 

For ii. for instance, assume that the sup of U on [x, y] is achieved at z = tx + (1 — t)y with 
t e]0, 1[. Then VU(z) ■ (y - x) = 0, so that 

1 <l-t) 2 



U(x)-U(z) = I (l-s)V 2 U(z + s(x-z))(x-z)-(x-z)ds^ inf a(Y) - — ■^- J -\x-y\ 
and similarly 



\Y-^\<R '1 



U(y) - U£±±) > U(y) - U(z) > inf a(Y) Ux - y' 2 
2 \ Y -z±y\<R 2 



Hence, for \x — y\ < 2R, 



U(z) - U(x) - U(y) < -U(?^l) - 2R 2 inf a(Y) < M. 

2 |y_£+2|< fl 

Proof of Proposition \2.2\ 

< Let (p be a strictly convex function on R n (with u = V such that L n Vip dfj, = J* K „ x d/i. 
First observe that 



\V<p(x) - V<p(y) -(x- y)\ z dfi(x) dfi(y) = 2 / \V<p(x) - x\ l dn(x) 
since, by assumption on (p, the difference is 

2 / (V(p(x) - x) dfi(x) 2 = 0. 

Then, by [H Lem. 5.1], 

(VW(x)-VW(y))-(x-y) > j\x - y\ 2 (9) 

if |x| ^ 2R or \y\ ^ 2i?. In view of this result we let 

X = {(x, y) £ R 2n ; \x - y\ < 2R, \V<p(x) - V<p(y)\ < 2R}. 
1. First of all, by convexity of W and ([9]), 

(VW(Vp(x) - Vp(y)) - Vlf(x - y)) • (Vy>(x) - V^(y) - (x - y)) dfi{x)dn{y) 



^ / (VW(V^(x) - V<p(y)) - VW(x - y)) ■ (V^(x) - V(f(y) - (x - y)) dn(x)dn(y) 

JR 2n \X 

> - / |Vy>(x) - V^(y) - (x - y)\ 2 dfi(x)dfi(y). 
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2. Then for all x and y, written as y = x + r 9 with r ^ and 9 G §> n , 
Vy?(y) - Vp(z) - (y - a;) = / [VV(a; + rte)--f]rfl(li 

JO 

We let H = VV(x + rt9) and write H- I = [H 1 / 2 - H* 1 / 2 ]!! 1 / 2 , so that 

< WH 1 / 2 - H'^W \H l ' 2 6\. 

Hence 

|Vp(j,) - Vip(x) -{y- x )\ 2 <r I' 1 \\H^ 2 - H^f e - u( - x+rte Ut C \ H 1 / 2 9\ 2 e u ^ x+rt ^ r dt 

Jo Jo 

by the Cauchy-Schwarz inequality. On the one hand, letting D = Aip + A</?*(V</?) — 2n, 

_ ff -i/ 2| |2 = gup ([g ~ 2/ + f" 1]3;) ' :E < irace(tf - 2/ + tf- 1 ) = D(x + rt9). 

x \X\ 

On the other hand 

f \H 1/2 9\ 2 rdt= [ V 2 p(x + rt9)(r9)-9dt = (Vip(y)-Vip(x))-9<2R 
Jo Jo 

if (x, y) G X. Hence 

|Vp(y) - Vv9(x) - (y-x)| 2 < AR 2 sup e* 7 ^ [ D(x + rt6)e- u( - x+rte Ut 

ze[x,y] Jo 

for all (x,y) G X, so that 



|V^(y) - Vip(x) - (y-x)\ dfj,(x)dfj,(y) 
<4R 2 [ e- u( - x) dx [ dy sup e^e"^ f 1 D(x + t(y - x)) e~ u(x+t( - y " x)) dt 

ill" J|j/-x|<2_R «£[a!,i/] JO 

<4R 2 e M / dz / f D{x + t{y - x)) e- u{x+ti ~ y ~ x)) dt 

il" J\y~x\<2R Jo 

by (jHJ). Now, for fixed f G [0, 1], the change of variables (x,y) h-» (u,it) = (x + t(y — x),y — x) 
has unit Jacobian, so this is equal to 



4i? 2 e M / 



dw / duD(v)e- u(v) = c [ D{v)d^{i 

J\u\<2R JR" 



for a constant c = c(-R, M, n) = 4R 2+n e M c n where c n is the volume of the unit ball in M. n . 
3. Collecting the terms in 1. and 2. concludes the proof with C = 2(j| +4c n R 2+n e M )~ 1 . \> 
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2.2 In presence of an exterior potential 

We saw in Section [T] how an exterior potential V can induce the convergence of all solutions to 
a unique equilibrium, and not only to the unique equilibrium with same center of mass as the 
initial datum of the solution to be considered. 

If W strictly convex (but at 0), and uniformly at infinity, and if V is strictly convex (but 
at 0), then polynomial convergence holds to a unique equilibrium fj,^ (see Section Q] and 
Th. 2.3]), and even exponential convergence, but with a rate depending on the free energy F 
of the initial datum (see Th. 2.5]). Following Theorem 12.11 (for V = 0), one may wonder 
whether this convergence is actually uniform in the initial datum, given by 

W 2 (jh, ^oo) < e- ct W 2 (fi ,^oc), t > 
for all solution (pt)t- But, according to Section [TJ this estimate is based on the inequality 

CWi(y,n) < Jv,w(v\ri ( 10 ) 

for the measure (jl = and all measures u; and this inequality does not hold if V is only 
assumed to be strictly convex. For instance: 

Lemma 2.4 Let fi € i^O^) an d V be a C 2 map on R, with V" bounded and V" — >+oo 0. Then 
there is no constant C > such that (|l(jp hold for all v. 

< We prove that fjlOf) does not hold for the translations v = <p'#n where <p'(x) = x + M 
where M —¥ +oo, that is, that there is no C > such that 

CM < f (V'(x + M) — V{x)) dfi(x) (11) 



for all M > 0. For that, we let R to be fixed later on, and bound the right-hand side in (jlip by 



/ 



-R 

\V'{x)\dp{x)+ / \V'(x + M)\dfi{x)+ I \V'(x + M)\dn(x). 

R 



First of all, since \V"\ < A, then |V'(a;)| < |V'(0)| + A\x\ so the first integral is finite 
(uniformly in M), and the second one is bounded by 

r—R r-R 

aM dfi(x)+ (\V'(0)\+A\x\)dfi(x). 



Now, for fixed e > 0, we take R such that this is bounded by (M + l)e for all M. Then we take 
Mq such that |V"(a;)| < e for x > M . For all M > M + R the third integral is bounded by 



\V'(M )\ + e(x + M - M )) d(i(x) < \V'(M )\ +e(m + M + I \x\ dfi(xj) . 



»+oo 
-R 

Collecting all terms we conclude that the full right-hand side in (|lip is < Ae for large M. > 
Lemma 12.41 only gives an instance of condition on V for (I10p not to hold. For example, the 



assumption V" bounded can be replaced by the doubling condition ([2]) for V' and J \ V'\ dfi < oo. 
Similarly, the assumption V" — > +OQ can be replaced by / + °° \ V"\ dfi < oo: in this case we use 

r+oo f+oo rx+M r+oo r 

I (V'(x + M) — V'{x)) d(j,(x) = / V"(t)dtdfi(x) < / \V"(t)\ dt dpi < C. 

-R J-R Jx J-R JR 
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Hence we can not expect a uniform rate of convergence to equilibrium for degenerately 
convex potentials. Our method is however able to recover an exponential convergence with a 
rate depending on the initial datum, as in [7J Th. 2.5]: 

Theorem 2.5 Assume that V is convex on W 1 with V 2 V(x) definite positive on \x\ ^ R and 
f e~ v dx < oo, and that W is convex with V 2 W(x) ^ K for \x\ ^ R. Then there exists a unique 
stationary solution /j,^ G (M") to (JTJ) . Moreover for all M there exists a positive constant C 
such that for all solutions (nt)t with J Rn \x\ 2 d/j,o(x) < M 

W 2 (nt, Moo) < e~ ct W 2 (no, Hoc), t ^ 0. 



< First, Proposition IA.ll ii. ensures the existence of a stationary measure /Uqo £ P2(M n ) 
which has a density satisfying /iioo(x) = Z~ 1 e~ v ( x ~ , ~ w *^° ( x ^ ; here Z is the normalizing constant. 
We just mention that the assumptions on V are satisfied by [2j Lem. 2.2] for instance. 

Then, by direct estimates on the propagation of the second moment, for all solutions (fJ-t)t 
with J Rn |x| 2 d^o(^) < M there is a constant N, depending only on V, W and M such that 

sup / \x\ 2 d[it{x) < N. 



Moreover, for v = Vy^/ioo with J Rn \x\ 2 dv(x) < N, we first write 



| V(p(x) — x\ 2 dfXoo = (V (p(x) — x) dfi 



+ ^ / 1 1^0*0 ~ ^fiv) ~{x- y)\ 2 dfi 00 (x)dfi 00 (y). 
By Proposition 12.61 below, applied with the constant ./V and the measure /Xqo, 

f(y<p(x) - x)dfi oa <\ I \V(p(x) - x\ 2 dfioo(x) + C J V)0 (^\fioo)- 
Then, by the proof of Proposition 12.21 there exists C\, depending only on V and W, such that 

\Vip(x) - Vtp(y) - (x - y)\ 2 dfx^x) d^oo(y) < C\J QjW {v\l*>oo)- 
Hence there exists a new positive constant C, depending only on V, W and M, such that 

CW 2 (fit, Moo) < Jv,w(lM\ Moo ) 

for all t. This proves the estimate on the convergence to /Xqo again by Proposition 11.11 and the 
Gronwall lemma. Uniqueness of the stationary solution in P 2 (W l ) follows. t> 

Proposition 2.6 Let V be a C 2 convex map on M n with V 2 V(x) definite positive on \x\ ^ R, 
and dfi(x) = e~ u<yX ^ dx be a probability measure on W 1 with U continuous. Then for all N there 
exists a constant C such that for all C 2 strictly convex map (p on W 1 with / R n|V99(x)| 2 dfj, < N 



V ip(x) dfj,(x) — / xdn(x 



2 1 
< - 
~ 2 



\V<p(x) - x\ 2 dfi(x) + C J Vfi (Vif#n\n). 
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<l Let S ^ 3R to be fixed later on. Since V is C 2 and V 2 V(x) is definite positive on the 
compact set R < \x\ < S, there exists K = K(S) > such that V 2 V(x) ^ K for all R < \x\ < S. 
Then, following [H Lem. 5.1], 

(VV(y) - VV(x)) -(y-x)> —\x - y\ 2 

if \x\ < S, \y\ < S and if |x| ^ 2R or \y\ ^ 2R; indeed one only need to take into account the 
values of V 2 V on the ball of radius S. 

Then we let ip be a given C 2 strictly convex map on M n and let 

X\ ={i£ R n , \x\ < S, \V<p(x)\ < S, \x\ ^ 2R or \<p(x)\ ^ 2R}. 
1. First of all, by convexity of V, the above remark and the Cauchy-Schwarz inequality, 
/ (W(V^) - W) • (Vp - x)d/i > [ (W(Vy>) - W) • {Vtp - x)dy. 

jR n JXi 

K 



K 



|V(/3 — x\ d/j, ^ 



(V(p — x)dn 



2. Then, on W 1 \ X u and letting X 2 = {x G R n , \x\ < 2R, \Vip(x)\ < 2R}, 



(Vip — x)d/i 



< 3 



\\7ip—x\d/j,) + 3 



x\^S 



X?ip{x)\^S 



\V(p— x\dfj,) +3( / \V(p — x\dfj, 



x 2 



By the Cauchy-Schwarz and Markov inequalities, the first term is bounded from above by 

J \V(p — x\ 2 dfi /j,[x, \x\ ^ S] < I I Vip — x\ 2 dfi — ^ J \x\ 2 dfi(x) 
J\x\^s Ju n Jm. n 

and the second one by 

/ \V(p-x\ 2 dfx^ f \V^(x)\ 2 dfi(x). 

Then, following the proof of [H Prop. 3.5], there exists a constant C, depending on V and 
/i only on the ball of radius 3R, such that 

/ |V(^ - x\ 2 dfi <C I (Ap + Aip*(Vp) -2n+ (W(Vy>) - W) ■ {Vtp - x)) dp. 

Jx 2 J\x\<3R 

This is in turn bounded by the corresponding integral on W 1 , which is Jv$(Vip#fi\fi), since 
both terms in the integrand are nonnegative. 

3. Collecting all terms we obtain 
6 



Vipd/j,— / xd/j, 

■n J^n 

Then we let S = max 

6 r 

S 2 



\x\ 2 dfi + N 



d»+ / \Vcp\ 2 di2+(—+6c)j v ,o(Vcp#n\n). 
} 



so that 



xfd(i+ / \V(f \ dfi 



1 

< - 
~ 2 



if J \V<p\ 2 dfi < N, concluding the proof with a C depending on V,fi and M through K(S). > 
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3 Non convex examples 

In this section we deal with potentials V and W for which the convergence rate to equilibrium 
is driven by V rather than by W. Our first result is more qualitative rather than quantitative. 

Theorem 3.1 Assume that V and W are C 2 convex maps and that there exist R > and K > 
such that for all \x\ > R, 

V 2 V(x) > K. 

Then there exists a unique stationary solution fj,^ £ P2(W n ) to flU, and a constant C such that 
for all solution (fJ.t)t of ([T]), 

Wadit.Atoo) < e -Cft W 2 (A«o,A»oo)i t> 0. 

In the first section (second example) we saw that only polynomial decay in contraction is 
known in this context, and only when the convexity degenerates at some points, for instance for 
V(x) = \x\ 4 . 

Proof 

< Existence of a stationary solution in P2(M. n ) which has a positive density satisfying 
Hoo = Z~ 1 e~ v ~ w * fJi °° is given by Proposition lA.il Hi, with any a < K and —a < b < 0. Then, 
by [H Prop. 3.5], there exists C > such that 

CWf (jh, Moo) < Jv,o{lh\ 

for all solution (fit)t- Moreover W is convex, so Jy.o < Jv,w- This proves the convergence bound 
by Proposition 11.11 Uniqueness of the stationary solution in P2(^ n ) follows. > 



Remark 3.2 The case of a double well potential for V is considered by J. Tugaut in \FT\ [73|/ . 
where the long time behavior is studied by a compactness argument, hence without rate. Let 
us now explain how Theorem lff.il extends to this case, for instance for V s {x) = x 4 — ex 2 and 
W{x) = |x| 3 in M. 

First of all, a stationary solution, solution of /i^, = e~^ E ~^ l '° a * w jZ e , exists by Proposi- 
tion \A.l\ Hi. Then one can then easily build a cut-off function tp such that V £ tp is C 2 , convex, 
satisfies (y e ip)" ^ K > outside a centered ball, uniformly in e £ [0,1], and is such that 
\\(V £ (1 — tp))"\\oo converges to as e — )• 0. Then, by Prop. 3.5], the measure //^ satisfies 
a WJys^ t Q inequality with a constant C > uniformly in e E [0, 1] (here we use that fWdula 
and Z £ are bounded uniformly in e). Now the perturbation proposition J^l Prop. 3.8] ensures 
that fj,^ satisfies a W Jv s o inequality, for sufficiently small e, hence a WJy s ,w inequality since 
W is convex. Here we say that a measure fj, satisfies a WJy,w inequality is the inequality (fTU|) 
holds for a positive constant C and all v . 

The smallness condition on e is necessary since, according to \18] . there exist several 
stationary solutions for large e. 

The following theorem provides the first examples of exponential convergence to equilibrium 
for the granular media equation, with both potentials non convex. 



11 



Theorem 3.3 Assume that 

• e~ v S P2(M. n ) and there exist a£l and C > such that V 2 V > a, and for all v 

Wi(v,e- v )<±J v , (v\e~ v ); (12) 

• there exist K ^ and /3 < such that sup < K and V 2 W > (3. 

Then there exists a unique stationary solution fi^ 6 P2(]R ra ) to §]]). Moreover, for all solution 
(m)t of and with C = (C - a)e~ 2K + a + /3, 

W2(fH,^oo) < e' dt W 2 (l^ 0l Hoo), t^O. 

Assumption (|12p on the measure e~ v has been studied in [1] under the name of WJ(C) 
inequality; there practical criteria have been given for the inequality to hold. Observe that we 
can always assume that C ^ a since, if a ^ 0, then \x satisfies a WJ(a) inequality. 

Proof 

< Existence of a stationary solution fj,^ in P2(M. n ) which has a positive density satisfying 
Hoo = Z~ l e is given by Proposition I A. 1\ iv; indeed, by [U Cor. 3.11], assumption (fT2l) 

on the measure e~ v implies the Talagrand inequality (I13p between the Wasserstein distance and 
the relative entropy (also called WH or T 2 , see [HH Chap. 22]), with the same C. 

Then we let dfi(x) = e~ v ^ dx and use the convexity assumptions on V and W, the bound 
on W and the sign conditions on /3 and C — a to get, for all v = V</?#Moo> 

Jv,w(v\lJ>oo) 

> e — J (Aip + A(p*{V<p) - 2n) dfi + J [(W(V<£>) - VV) ■ {Vip - x) - a\Vip - x\ 2 } d^ 
+a j \V<p- x\ 2 dHoo + ^ ^ \Vf(x) -Vip(y) - (x - y)\ 2 dfi 00 (x)dfj, 00 (y) 

—K r —K r 

> —I (A^ + A(p* t (V^) -2n)dfi+ ^- / [{VV(V<p) - VV)-(V(f - x) - a\V<p - x\ 2 } dpt 



+(a + /3)J |V^-x| 2 d/ioo -p j\v^-x)d^L 00 2 

> (C-a)^r- J \V<p-x\ 2 dfx + {a + 13) J \Vip-x\ 2 dfi 

> C J \Vip - x\ 2 dfiocix) = CWi(u, Moo ). 



> 



A Existence of stationary solutions 

The existence of a minimizer of F has been proved by R. J. McCann [16] for strictly convex 
or radially symmetric convex interaction potentials W (and V = 0). We adapt his classical 
compactness-lower semicontinuity argument to our diverse cases: 

Proposition A.l The map F : P2(W l ) ->1U {+00} defined by (|3j) for absolutely continuous 
measures and by +00 otherwise achieve its minimum in each of the following cases: 
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i. V = 0, W is convex and W(x) ^ b\x\ 2 — b' for 6, b' > 0; 

ii. V(x) ^ a\x\ — a' and W{x) ^ b\x\ 2 — b' for a, a', b, b' > 0; 

Hi. V(x) ^ a|x| 2 — a' and W(x) ^ 6|x| 2 — b' for b' , a, a' > 0, b > —a; 

iv. W is bounded from below and e~ v £ i-^O^™) satisfies a Talagrand transportation inequality 

2 



W^u,e~ v )< 



C 



I 



v log v dx + Vdu), ueP 2 



(13) 



Then, as in [7], a minimizer of F has a positive density on W 1 satisfying 

log floo + V + W * floe = A Gl. 

Proo/ 

<l First of all, mfp 2 ( R n\ F < +oo since F(fi) < +oo for /x the Lebesgue measure on [0,1]™ 
for instance. Let then {[i p ) p € P2(M. n ) be a minimizing sequence, and assume for a while that 
J \x\ 2 d[i p is bounded. Then (/i p ) p is tight, so up to a subsequence admits a limit for the 
narrow convergence by the Prohorov Theorem. Moreover J \x\ 2 dfioo < liminf p J \x\ 2 d^L p < +oo 
so /Xqo G P2(^ n )- Finally /Uqo minimizes -F on P2(IR n ) by lower semicontinuity. 

It remains now to bound J \x\ 2 d/j, p by F(/jl p ) in each case: 

For i., as in |16J . let V<y? p transport /U p onto ^ p (— .) and let /i p = #/x p for I the identity 

map. Now IF is convex, so F is displacement convex, so that F{fl p ) < (F(/i p ) + F(n p (— .)))/2 = 
F(fi p ) and (/i p ) is also a minimizing sequence. Moreover J xdfi p = so 



\x\ 2 dfi p 



\x - y\ 2 dfi p (x)dfj, p (y) < 



-(W(x -y) + b')dfiJx)dfip(y) < 



+ 



2b 



For ii. we observe that 







2- 


y I 


j \x\ 2 dfi p — 


J xdfi p 







b< 



hence f \x\d/j, p is bounded by the second inequality, and then J \x\ 2 dfj, p by the first one. 
For Hi. we similarly observe, and by discussing on the sign of b, that 

F(n p ) ^ a f \x\ 2 dfi p — a'+— f I ' (b\x—y\ 2 —b')dfi p (x)dfj, p (y) ^ (a+min(6, 0)) / \x\ 2 dfi p — a— — -■ 



For iv. we notice that 



J \x\ 2 d^ p < 2Wi(fi p ,e- v ) + 2 j 



2„-V 



x\ e 



< 



C 
4 



fj, p log n P dx + j V dfi p J + 2 
1 



I |2 -V 

\x\ e 



inf W + 2 



|2 -V 

x\ e . 



> 
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